\(\int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx\) [348]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [B] (verified)
   Fricas [B] (verification not implemented)
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 25, antiderivative size = 91 \[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\frac {(a+2 b) \text {arctanh}\left (\frac {\sqrt {a+b} \sin (e+f x)}{\sqrt {a+b \sin ^2(e+f x)}}\right )}{2 (a+b)^{3/2} f}+\frac {\sec (e+f x) \sqrt {a+b \sin ^2(e+f x)} \tan (e+f x)}{2 (a+b) f} \]

[Out]

1/2*(a+2*b)*arctanh(sin(f*x+e)*(a+b)^(1/2)/(a+b*sin(f*x+e)^2)^(1/2))/(a+b)^(3/2)/f+1/2*sec(f*x+e)*(a+b*sin(f*x
+e)^2)^(1/2)*tan(f*x+e)/(a+b)/f

Rubi [A] (verified)

Time = 0.13 (sec) , antiderivative size = 91, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.160, Rules used = {3269, 390, 385, 212} \[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\frac {(a+2 b) \text {arctanh}\left (\frac {\sqrt {a+b} \sin (e+f x)}{\sqrt {a+b \sin ^2(e+f x)}}\right )}{2 f (a+b)^{3/2}}+\frac {\tan (e+f x) \sec (e+f x) \sqrt {a+b \sin ^2(e+f x)}}{2 f (a+b)} \]

[In]

Int[Sec[e + f*x]^3/Sqrt[a + b*Sin[e + f*x]^2],x]

[Out]

((a + 2*b)*ArcTanh[(Sqrt[a + b]*Sin[e + f*x])/Sqrt[a + b*Sin[e + f*x]^2]])/(2*(a + b)^(3/2)*f) + (Sec[e + f*x]
*Sqrt[a + b*Sin[e + f*x]^2]*Tan[e + f*x])/(2*(a + b)*f)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 385

Int[((a_) + (b_.)*(x_)^(n_))^(p_)/((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Subst[Int[1/(c - (b*c - a*d)*x^n), x]
, x, x/(a + b*x^n)^(1/n)] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - a*d, 0] && EqQ[n*p + 1, 0] && IntegerQ[n]

Rule 390

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[(-b)*x*(a + b*x^n)^(p + 1)*
((c + d*x^n)^(q + 1)/(a*n*(p + 1)*(b*c - a*d))), x] + Dist[(b*c + n*(p + 1)*(b*c - a*d))/(a*n*(p + 1)*(b*c - a
*d)), Int[(a + b*x^n)^(p + 1)*(c + d*x^n)^q, x], x] /; FreeQ[{a, b, c, d, n, q}, x] && NeQ[b*c - a*d, 0] && Eq
Q[n*(p + q + 2) + 1, 0] && (LtQ[p, -1] ||  !LtQ[q, -1]) && NeQ[p, -1]

Rule 3269

Int[cos[(e_.) + (f_.)*(x_)]^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(p_.), x_Symbol] :> With[{ff = Free
Factors[Sin[e + f*x], x]}, Dist[ff/f, Subst[Int[(1 - ff^2*x^2)^((m - 1)/2)*(a + b*ff^2*x^2)^p, x], x, Sin[e +
f*x]/ff], x]] /; FreeQ[{a, b, e, f, p}, x] && IntegerQ[(m - 1)/2]

Rubi steps \begin{align*} \text {integral}& = \frac {\text {Subst}\left (\int \frac {1}{\left (1-x^2\right )^2 \sqrt {a+b x^2}} \, dx,x,\sin (e+f x)\right )}{f} \\ & = \frac {\sec (e+f x) \sqrt {a+b \sin ^2(e+f x)} \tan (e+f x)}{2 (a+b) f}+\frac {(a+2 b) \text {Subst}\left (\int \frac {1}{\left (1-x^2\right ) \sqrt {a+b x^2}} \, dx,x,\sin (e+f x)\right )}{2 (a+b) f} \\ & = \frac {\sec (e+f x) \sqrt {a+b \sin ^2(e+f x)} \tan (e+f x)}{2 (a+b) f}+\frac {(a+2 b) \text {Subst}\left (\int \frac {1}{1-(a+b) x^2} \, dx,x,\frac {\sin (e+f x)}{\sqrt {a+b \sin ^2(e+f x)}}\right )}{2 (a+b) f} \\ & = \frac {(a+2 b) \text {arctanh}\left (\frac {\sqrt {a+b} \sin (e+f x)}{\sqrt {a+b \sin ^2(e+f x)}}\right )}{2 (a+b)^{3/2} f}+\frac {\sec (e+f x) \sqrt {a+b \sin ^2(e+f x)} \tan (e+f x)}{2 (a+b) f} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 5 vs. order 3 in optimal.

Time = 10.14 (sec) , antiderivative size = 408, normalized size of antiderivative = 4.48 \[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\frac {\sec ^3(e+f x) \left (1+\frac {b \sin ^2(e+f x)}{a}\right ) \tan (e+f x) \left (45 a \arcsin \left (\sqrt {-\frac {(a+b) \tan ^2(e+f x)}{a}}\right )+30 b \arcsin \left (\sqrt {-\frac {(a+b) \tan ^2(e+f x)}{a}}\right ) \sin ^2(e+f x)+16 a \operatorname {Hypergeometric2F1}\left (2,3,\frac {7}{2},-\frac {(a+b) \tan ^2(e+f x)}{a}\right ) \sqrt {\frac {\sec ^2(e+f x) \left (a+b \sin ^2(e+f x)\right )}{a}} \left (-\frac {(a+b) \tan ^2(e+f x)}{a}\right )^{5/2}+16 b \operatorname {Hypergeometric2F1}\left (2,3,\frac {7}{2},-\frac {(a+b) \tan ^2(e+f x)}{a}\right ) \sin ^2(e+f x) \sqrt {\frac {\sec ^2(e+f x) \left (a+b \sin ^2(e+f x)\right )}{a}} \left (-\frac {(a+b) \tan ^2(e+f x)}{a}\right )^{5/2}-45 a \sqrt {-\frac {(a+b) \sec ^2(e+f x) \left (a+b \sin ^2(e+f x)\right ) \tan ^2(e+f x)}{a^2}}-30 b \sin ^2(e+f x) \sqrt {-\frac {(a+b) \sec ^2(e+f x) \left (a+b \sin ^2(e+f x)\right ) \tan ^2(e+f x)}{a^2}}\right )}{30 a f \sqrt {a+b \sin ^2(e+f x)} \sqrt {\frac {\sec ^2(e+f x) \left (a+b \sin ^2(e+f x)\right )}{a}} \left (-\frac {(a+b) \tan ^2(e+f x)}{a}\right )^{3/2}} \]

[In]

Integrate[Sec[e + f*x]^3/Sqrt[a + b*Sin[e + f*x]^2],x]

[Out]

(Sec[e + f*x]^3*(1 + (b*Sin[e + f*x]^2)/a)*Tan[e + f*x]*(45*a*ArcSin[Sqrt[-(((a + b)*Tan[e + f*x]^2)/a)]] + 30
*b*ArcSin[Sqrt[-(((a + b)*Tan[e + f*x]^2)/a)]]*Sin[e + f*x]^2 + 16*a*Hypergeometric2F1[2, 3, 7/2, -(((a + b)*T
an[e + f*x]^2)/a)]*Sqrt[(Sec[e + f*x]^2*(a + b*Sin[e + f*x]^2))/a]*(-(((a + b)*Tan[e + f*x]^2)/a))^(5/2) + 16*
b*Hypergeometric2F1[2, 3, 7/2, -(((a + b)*Tan[e + f*x]^2)/a)]*Sin[e + f*x]^2*Sqrt[(Sec[e + f*x]^2*(a + b*Sin[e
 + f*x]^2))/a]*(-(((a + b)*Tan[e + f*x]^2)/a))^(5/2) - 45*a*Sqrt[-(((a + b)*Sec[e + f*x]^2*(a + b*Sin[e + f*x]
^2)*Tan[e + f*x]^2)/a^2)] - 30*b*Sin[e + f*x]^2*Sqrt[-(((a + b)*Sec[e + f*x]^2*(a + b*Sin[e + f*x]^2)*Tan[e +
f*x]^2)/a^2)]))/(30*a*f*Sqrt[a + b*Sin[e + f*x]^2]*Sqrt[(Sec[e + f*x]^2*(a + b*Sin[e + f*x]^2))/a]*(-(((a + b)
*Tan[e + f*x]^2)/a))^(3/2))

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(359\) vs. \(2(79)=158\).

Time = 1.05 (sec) , antiderivative size = 360, normalized size of antiderivative = 3.96

method result size
default \(\frac {-\left (\ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}-2 b \sin \left (f x +e \right )+2 a}{1+\sin \left (f x +e \right )}\right ) a^{2}+3 \ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}-2 b \sin \left (f x +e \right )+2 a}{1+\sin \left (f x +e \right )}\right ) a b +2 \ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}-2 b \sin \left (f x +e \right )+2 a}{1+\sin \left (f x +e \right )}\right ) b^{2}-\ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}+2 b \sin \left (f x +e \right )+2 a}{\sin \left (f x +e \right )-1}\right ) a^{2}-3 \ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}+2 b \sin \left (f x +e \right )+2 a}{\sin \left (f x +e \right )-1}\right ) a b -2 \ln \left (\frac {2 \sqrt {a +b}\, \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}+2 b \sin \left (f x +e \right )+2 a}{\sin \left (f x +e \right )-1}\right ) b^{2}\right ) \left (\cos ^{2}\left (f x +e \right )\right )+2 \sin \left (f x +e \right ) \sqrt {a +b -b \left (\cos ^{2}\left (f x +e \right )\right )}\, \left (a +b \right )^{\frac {3}{2}}}{4 \left (a +b \right )^{\frac {5}{2}} \cos \left (f x +e \right )^{2} f}\) \(360\)

[In]

int(sec(f*x+e)^3/(a+b*sin(f*x+e)^2)^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/4*(-(ln(2/(1+sin(f*x+e))*((a+b)^(1/2)*(a+b-b*cos(f*x+e)^2)^(1/2)-b*sin(f*x+e)+a))*a^2+3*ln(2/(1+sin(f*x+e))*
((a+b)^(1/2)*(a+b-b*cos(f*x+e)^2)^(1/2)-b*sin(f*x+e)+a))*a*b+2*ln(2/(1+sin(f*x+e))*((a+b)^(1/2)*(a+b-b*cos(f*x
+e)^2)^(1/2)-b*sin(f*x+e)+a))*b^2-ln(2/(sin(f*x+e)-1)*((a+b)^(1/2)*(a+b-b*cos(f*x+e)^2)^(1/2)+b*sin(f*x+e)+a))
*a^2-3*ln(2/(sin(f*x+e)-1)*((a+b)^(1/2)*(a+b-b*cos(f*x+e)^2)^(1/2)+b*sin(f*x+e)+a))*a*b-2*ln(2/(sin(f*x+e)-1)*
((a+b)^(1/2)*(a+b-b*cos(f*x+e)^2)^(1/2)+b*sin(f*x+e)+a))*b^2)*cos(f*x+e)^2+2*sin(f*x+e)*(a+b-b*cos(f*x+e)^2)^(
1/2)*(a+b)^(3/2))/(a+b)^(5/2)/cos(f*x+e)^2/f

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 167 vs. \(2 (79) = 158\).

Time = 0.43 (sec) , antiderivative size = 361, normalized size of antiderivative = 3.97 \[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\left [\frac {{\left (a + 2 \, b\right )} \sqrt {a + b} \cos \left (f x + e\right )^{2} \log \left (\frac {{\left (a^{2} + 8 \, a b + 8 \, b^{2}\right )} \cos \left (f x + e\right )^{4} - 8 \, {\left (a^{2} + 3 \, a b + 2 \, b^{2}\right )} \cos \left (f x + e\right )^{2} - 4 \, {\left ({\left (a + 2 \, b\right )} \cos \left (f x + e\right )^{2} - 2 \, a - 2 \, b\right )} \sqrt {-b \cos \left (f x + e\right )^{2} + a + b} \sqrt {a + b} \sin \left (f x + e\right ) + 8 \, a^{2} + 16 \, a b + 8 \, b^{2}}{\cos \left (f x + e\right )^{4}}\right ) + 4 \, \sqrt {-b \cos \left (f x + e\right )^{2} + a + b} {\left (a + b\right )} \sin \left (f x + e\right )}{8 \, {\left (a^{2} + 2 \, a b + b^{2}\right )} f \cos \left (f x + e\right )^{2}}, -\frac {{\left (a + 2 \, b\right )} \sqrt {-a - b} \arctan \left (\frac {{\left ({\left (a + 2 \, b\right )} \cos \left (f x + e\right )^{2} - 2 \, a - 2 \, b\right )} \sqrt {-b \cos \left (f x + e\right )^{2} + a + b} \sqrt {-a - b}}{2 \, {\left ({\left (a b + b^{2}\right )} \cos \left (f x + e\right )^{2} - a^{2} - 2 \, a b - b^{2}\right )} \sin \left (f x + e\right )}\right ) \cos \left (f x + e\right )^{2} - 2 \, \sqrt {-b \cos \left (f x + e\right )^{2} + a + b} {\left (a + b\right )} \sin \left (f x + e\right )}{4 \, {\left (a^{2} + 2 \, a b + b^{2}\right )} f \cos \left (f x + e\right )^{2}}\right ] \]

[In]

integrate(sec(f*x+e)^3/(a+b*sin(f*x+e)^2)^(1/2),x, algorithm="fricas")

[Out]

[1/8*((a + 2*b)*sqrt(a + b)*cos(f*x + e)^2*log(((a^2 + 8*a*b + 8*b^2)*cos(f*x + e)^4 - 8*(a^2 + 3*a*b + 2*b^2)
*cos(f*x + e)^2 - 4*((a + 2*b)*cos(f*x + e)^2 - 2*a - 2*b)*sqrt(-b*cos(f*x + e)^2 + a + b)*sqrt(a + b)*sin(f*x
 + e) + 8*a^2 + 16*a*b + 8*b^2)/cos(f*x + e)^4) + 4*sqrt(-b*cos(f*x + e)^2 + a + b)*(a + b)*sin(f*x + e))/((a^
2 + 2*a*b + b^2)*f*cos(f*x + e)^2), -1/4*((a + 2*b)*sqrt(-a - b)*arctan(1/2*((a + 2*b)*cos(f*x + e)^2 - 2*a -
2*b)*sqrt(-b*cos(f*x + e)^2 + a + b)*sqrt(-a - b)/(((a*b + b^2)*cos(f*x + e)^2 - a^2 - 2*a*b - b^2)*sin(f*x +
e)))*cos(f*x + e)^2 - 2*sqrt(-b*cos(f*x + e)^2 + a + b)*(a + b)*sin(f*x + e))/((a^2 + 2*a*b + b^2)*f*cos(f*x +
 e)^2)]

Sympy [F]

\[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\int \frac {\sec ^{3}{\left (e + f x \right )}}{\sqrt {a + b \sin ^{2}{\left (e + f x \right )}}}\, dx \]

[In]

integrate(sec(f*x+e)**3/(a+b*sin(f*x+e)**2)**(1/2),x)

[Out]

Integral(sec(e + f*x)**3/sqrt(a + b*sin(e + f*x)**2), x)

Maxima [F]

\[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\int { \frac {\sec \left (f x + e\right )^{3}}{\sqrt {b \sin \left (f x + e\right )^{2} + a}} \,d x } \]

[In]

integrate(sec(f*x+e)^3/(a+b*sin(f*x+e)^2)^(1/2),x, algorithm="maxima")

[Out]

integrate(sec(f*x + e)^3/sqrt(b*sin(f*x + e)^2 + a), x)

Giac [F]

\[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\int { \frac {\sec \left (f x + e\right )^{3}}{\sqrt {b \sin \left (f x + e\right )^{2} + a}} \,d x } \]

[In]

integrate(sec(f*x+e)^3/(a+b*sin(f*x+e)^2)^(1/2),x, algorithm="giac")

[Out]

sage0*x

Mupad [F(-1)]

Timed out. \[ \int \frac {\sec ^3(e+f x)}{\sqrt {a+b \sin ^2(e+f x)}} \, dx=\int \frac {1}{{\cos \left (e+f\,x\right )}^3\,\sqrt {b\,{\sin \left (e+f\,x\right )}^2+a}} \,d x \]

[In]

int(1/(cos(e + f*x)^3*(a + b*sin(e + f*x)^2)^(1/2)),x)

[Out]

int(1/(cos(e + f*x)^3*(a + b*sin(e + f*x)^2)^(1/2)), x)